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Abstract. We establish quantum dynamical lower bounds for a number of 
discrete one-dimensional Schrodinger operators. These dynamical bounds are 
derived from power-law upper bounds on the norms of transfer matrices. We 
develop further the approach from part I and study many examples. Particular 
focus is put on models with finitely or at most countably many exceptional 
energies for which one can prove power-law bounds on transfer matrices. The 
models discussed in this paper include substitution models, Sturmian mod- 
els, a hierarchical model, the prime model, and a class of moderately sparse 
potentials. 

1. Introduction 
Consider a discrete one-dimensional Schrodinger operator 

(1) Hi/jin) = V(" - 1) + + 1) + V{n)ip{n), 

in ^^(Z) or £^(N) (with a Dirichlet boundary condition). We are interested in 
proving lower bounds on the spreading of an initially localized wavepacket under 
the dynamics governed by H. That is, if we consider the initial state V', we ask how 
fast ilj{t) = exp(— spreads out. One is normally interested in initial states 
that are well localized. In the present paper we shall limit our attention to the case 

A typical quantity that is considered to measure the spreading of ij.'{t) is the 
following: Define 

(2) (|AQ(r) = ^|nraKT), 
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where 

(3) a{n,T)^-J^ e~''^^\{Sn,m)\' dt. 

Clearly, the faster {\X\^){T) grows, the faster ij;{t) spreads out, at least averaged 
in time. One typically wants to prove power-law lower bounds on and 
hence it is natural to define the following quantity: For p > 0, define the lower 
growth exponent (3^ (p) by 



log(|^iP(T) 



T^+oo logT 

There are presently two distinct approaches to proving lower bounds for (p) . 
The first goes back to works of Guarneri Combes and Last |^ and is 
based on a study of the Hausdorff dimension of the spectral measure /z^ associated 
with the pair (iJ, ip). Namely, we have the following bound: 

(4) >p•dimH(Aiv)• 
The Jitomirskaya-Last extension |15[ [TB| of Gilbert-Pearson theory allows for 
a convenient way of investigating dim//(/Lt^) and hence this approach has enjoyed 
some popularity (see, e.g., 01231 EHl for applications). 

On the other hand, this bound clearly gives nothing in the case of a zero- 
dimensional spectral measure, for example, in the case of a pure point measure. 
There are a number of models where one expects (or can prove) pure point spec- 
trum with strictly positive values for P^{p)- An example is given by the random 
dimer model; studied, for example, in El Ej ■ It is therefore desirable to have 
a way of proving lower bounds on the transport exponents which works for such 
models and, of course, whose input is easy to verify in concrete cases. Such an 
approach was developed in [5| (and employed in ^7] to prove the conjectured dy- 
namical lower bound for the random dimer model), and the present article is a 
continuation of that paper. The necessary input are power-law upper bounds on 
transfer matrices for certain energies. It may come as a surprise that dynamical 
bounds can be obtained if there is only one energy where one can exhibit a power- 
law bound for the transfer matrix. This is indeed necessary for models such as the 
random dimer model and related ones 7], where there are only a finite number of 
such energies. 

Another advantage of the approach from |H] over the bound 10} is the stability 
of its input with respect to perturbations of the potential V. It was noted in jHj 
that if its approach can be applied to a given model, then it can also be applied to 
all finitely supported perturbations of the given potential — and it gives the same 
dynamical bounds for the perturbed models. Such a stability is not true in general 
for bounds derived using Q). For example, it may happen that the addition of a 
finitely supported perturbation turns a given singular continuous spectral measure 
into a pure point measure; see [5] for many examples illustrating this phenomenon. 

In the general criterion was applied to three prominent models from one- 
dimensional quasicrystal theory, namely, the Fibonacci model, the period doubling 
model, and the Thue-Morse model. All these models can be generated by a substi- 
tution process. This allows one to study the growth of transfer matrix norms with 
the help of an associated dynamical system — the trace map — and this provides in 
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particular a very convenient way of verifying the input to the general dynamical 
criterion. 

In the present paper we will prove a more general version of the dynamical result 
from involving also the weight assigned by the spectral measure to the set of 
energies with power-law bounded transfer matrices. This gives stronger dynamical 
results in cases where such bounds hold for all energies in the spectrum, for example, 
models with Sturmian potentials. We shall also prove a stronger stability result. 
Namely, we will show that, for a fixed energy, the power-law bound is stable with 
respect to power-decaying perturbations. Here, the power-decay of the perturbation 
that we can allow depends on the transfer matrix power-law bound we start out 
with. Finally, we shall study a large number of examples and derive dynamical 
results for them by applying our main theorem, Theorem ^ below. The examples 
discussed in this paper include, in particular, generalizations of each of the three 
prominent substitution models studied in 0. 

The organization of the paper is as follows. In Section [5] we prove our main 
theorem which derives quantum dynamical lower bounds from power-law bounds 
on transfer matrices. Section |3| discusses the stability of such power-law bounds 
on transfer matrices with respect to power-decaying perturbations of the potential. 
Section 0] deals with a class of models that are "sparse" in a certain sense and 
which includes a variety of substitution models (in particular, generalizations of 
Fibonacci, period doubling, and Thue-Morse), the prime Schrodinger operator, and 
moderately sparse models which were studied by Zlatos in 35 . The hierarchical 
model, which was studied in detail by Kunz et al. in |23j from a spectral point of 
view, will then be considered in Section [3 Finally, we present results for Sturmian 
models (studied, e.g., in QJISIE]; see also the reviews 01123) in Sectional 

2. A Quantum Dynamical Lower Bound Derived From Power-Law 



In this section, we prove a more general version of the main result from 8 . The 
general idea of proof is the same and the result derives lower bounds on the dy- 
namical quantity (p) from power-law bounds on transfer matrices. However, the 
result established in this section gives improved bounds in many cases, in particular 
in the case of Sturmian potentials discussed later in the paper. 

Recall the notion of a transfer matrix. Consider for some i? G M, a solution (j) of 
the difference equation 



Denote $(n) = {4>{n -|- 1) , (j){n))'^ . The transfer matrix T{n,m; E) is defined by 
requiring 

$(n) = T{n,m;E)<^{m) 
for every solution (/) of (0). It is straightforward to verify that for n > m, 
T{n, m; E) = T(y{n)]E) x • • • x T{V{m +iy,E), 



Transfer Matrix Bounds 



(5) 



(j){n + 1) + (j){n 



1) + y(n)(/)(n) = E(f>{Ti). 



where 




and similarly for n < m. 

With this notation at hand we can now state: 
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Theorem 1. The following statements hold: 

(a) Suppose that for some K > Q, C > Q, a > the following condition holds: For 
any N > large enough, there exists a nonempty Borel set A{N) c M such that 

A{N) C [-K,K] and 

(6) ||r(n, to; E)\\ < CN"" \fE G A{N), V n, to : |n| < TV, |to| < N 

{resp., with I < n < N, I < m < N m the case off{N)). Let N(T) = T^/(^+») 
and let, for j = 1, 2, Bj{T) be the j / T -neighborhood of the set A{N{T)): 

Bj{T) = e K : e A{N{T)), \E - E'\ < j/T}. 

Denote by F{z) the Borel transform of the spectral measure of the state ip = 5i: 

d^i{x) 



F{E + is 



ie)= [ 
Jr 



X- {E + ie)' 



Then for the initial state tjj = 6i and allT > 1 large enough, the following bound 
holds: 

(7) P{T) = J2 ^) ^ ^N'~^"{T) [ dE{l + Wf{E + is)). 

n:\n\>N{T) ^ •' B2{T) 

In particular, 

(8) P(T) > ^iVi-2"(T)(|i3i(T)| + 

where \B\ denotes the Lebesgue measure. This gives the following bound for the 
time- averaged moments: 

(9) > ^NP+'-^'^{T)i\B,{T)\+^^{B,{T))). 

(b) Suppose that there exists a set A c [—K,K] of positive measure n{A) > such 
that 

\\T{n,m;E)\\<C{\nr + \mn 
for all E G A, n, to. Then 

(c) Assume that 

\\T{n,m;Eo)\\<C{Eo){\n\" + \mn 
for some Eq, uniformly in n,m, then 

(11) {\X\l)iT) > CT"^ +fx{[Eo- T-\ Eo + T-'])) . 

Assume moreover that Eq is an eigenvalue {possible only if a > 1/2) , .so that there 
exists ifj € 7^ such that Hip = Eotjj- Suppose that V'(l) 7^ {this is always 

true in the case of£^{N)). Then 
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Proof. As in |H] we shall consider the case of £^(Z), because for ^^(N), the proof is 
similar but simpler. The main part of the proof is virtually identical with that of 
[5]. For the sake of completeness we shall briefly recall the main lines. 
The starting point is the Parseval equality: 

a(n,T) = - / e~^''^\{5n,eM-itH)5i)\''dt = — / \{5n, R{E + ie)6i)\'' dE, 
^ Jo Jm 

where R{z) = {H~zl)~^ and e = l/T. For z — E+ie, e > 0, we define (p = R{z)Si, 
^{n) — {4>{n + 1), (t){n)Y' . For each n > 1, one has the inequality 

(13) ||$(n)||>||T(n,l;z)|ri||$(l)|| 
and for each n < 0, 

(14) ||$(n)||>||T(n,0;z)|ri||<i>(0)||. 

An upper bound for the norm of the transfer matrix with complex z is obtained 
using condition JH)) and Lemma 2.1]. Namely, let us fix some T>1, e ^ l/T 
and define N = NiT) = T^/(.'^+°'). Then for every E e B2{T) and I < n < N 

(15) \\T{n,l;E + ie)\\<DN°', 

where D — Cexp(3C), and C is the constant from ®. A similar bound holds 
for negative values of n. Using the bounds H13|l - H15|) . one shows that for every 
E £ B2{T), 

(16) \{Sn,R{E+te)s,)\^>cN^-''^{\m\^ + m)\' + \m\') 

n:\n\>N/2 

with uniform constant c > 0. It was shown in |Sj that under the conditions of the 
theorem one always has 

1^(0)1 + 1,^(1)1 + |<^(2)| > c> 

with uniform constant. What one can also observe (and this is a new point) is the 
fact that 

0(1) - {R{z)Si,Si) ^F{z), 
where F{z) is the Borel transform of the spectral measure corresponding to the 
pair {H,6i). Therefore, it follows from H16|l that 

^ \{6n,RiE + ie)6i)\^ > cN^-^'^{l + WF{E + ie)). 

n:\n\>N/2 

Integrating this bound over E e B2{T), one proves ((T)). Next, one observes that 
l + ImV(z) > 2ImF(z). For any set S, denote by the e-neighborhood of S. 
Following T^, one can see that 

edE 



/ lmF{E + ie)dE ^ / dfj.{x) j 
Js, Jr Js 

> / dfi{x) 



5. i^-Ey 

edu 



Taking S = Bi{T), we prove 0. The bound immediately follows. 

To prove part (b), one just takes A{N) — A for every N. Since ^(i?i(T)) > 
H{A{N{T)) = fi{A) > 0, the result follows from the bound ©. 
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The bound (fTT|l of part (c) follows directly from taking A{N) = {Eq} for 
every N. Finally, to prove the second part of (c), we go back to Q to obtain 

{\X\l)iT) > ^N^+^-'^iT) [ WF{E + ^e)dE, 

where B2{T) = [E(^~2£,E(^ + 2£\. Under the condition ipil) ^ 0, one has /i({£'o}) > 
0. Thus, 

'^^^""^''^^jE-irTe^- 
Integration over B2 {T) yields . □ 

Remark. Part (b) of Theorem ^ remains true if 

(17) \\T{n,m-E)\\<C{E){\nr + \mn 

for all n, m and E G A with C{E) < 00 for /x-almost every E. To prove this, it is 
sufficient to take a smaller set A' C A oi positive measure where C{E) < C < 00. 
The bound (|10|l should be compared with the well-known result of jl5[ : If H17|) 
holds for some a G [0, 1/2) on a set A of positive /^-measure, then the restriction of 
^ to A is 1 — 2a-continuous. In particular, 

/35;(p)>p(l-2a). 

This bound is better than IjlUfl for small p, but for p large enough, (|10|l is always 
better. Moreover, 1)10(1 holds also if a > 1/2. 

3. Stability With Respect to Power-Decaying Perturbations 

In this section, we discuss the stability of the crucial input to our dynamical 
bounds, power-law bounds on transfer matrices, with respect to perturbations of 
the potential. It is easy to see, and was noted in O Corollary 1.3], that finitely 
supported perturbations of the potential cannot destroy such a power-law bound. 
Here we strengthen this to stability with respect to power-decaying perturbations, 
where the allowed power depends on the bound we can prove for the unperturbed 
problem. 

Theorem 2. Assume that for some energy E and some constant Ci, the transfer 
matrices T associated with IS. -\-V satisfy 

(18) \\T{n, m; £')|| < Ci\n~ m|" for every n,m G Z with nm > 0. 
As.sume further that, for some e > 0, the perturbation W satisfies 

(19) \Win)\ < 02(1 + |n|)~i-2"-^ for every n e Z. 
Then the transfer matrices T' associated with A + V + W satisfy 

(20) \\T'{n, m; E)\\ < Csjn — m|" for every ri, m G Z with nm > 0. 

Proof. We present the proof in the special case where we assume (II8II only for n > 
and m — and then prove (|20|l for n > and m = 0. A slight variation of the 
argument below works for general ri, m G Z with nm > (with a uniform constant 
C3 in 

Our strategy will be to work with solutions and employ a general perturbation 
method developed by Kiselev et al. in '21'. 

Consider the unperturbed equation (jSJ and the perturbed equation 

(21) tP{n + ipin - 1) + [V{n) + W{n)]%l;{n) = E%l;{n). 



QUANTUM DYNAMICS IN ONE DIMENSION II. 



7 



Note that the transfer matrix T'{n, 0; E) is given by 

iJj-Din + l) i>-N{n- 



V'N(n) 



= /. 



T'{n,0;E) = 

where ipu , V'n solve lPT|) and obey 

f Vd(1) ^n(I) 
^ V-dW i^n(0) 

Fix a complex reference solution (/> of lOl- For example, we could set = +«0n, 
where (/)d, ^n solve (jSJ and have the same initial conditions as "^d, V'n- By (|18|l we 
have 

(22) |(/.(n)| < qnr. 

Let t/; be one of the basic solutions ^d, ^n of (|21|l . Define p{n) by 



■^(n — 1) 



2i 
Im 



p(n) 



p{n) 



{n-1) 

(f){n) 
An -I) 



-pin) 



1) 



Write 



2T4^(n) 



(n) and in polar coordinates, 

cj){n) = |0(n)|e''^(") pin) 

and define 

6* (n) = 77(71) + 7(71) and U{n) 

UJ 

where iuj is the Wronskian of cf) and (/), that is, 

2ilm(0(7i + l)0(n)) = zcj for every n. 

Clearly, the assertion of the theorem follows if we can show that Rin) remains 
bounded as |?t.| — > 00. The key identity (equation (45) in |2J) is the following: 

(23) Rin + 1)2 = i?(n)2[l + Uin) sin(26l(7i)) + [/(ri)^ sin^ (61(71))]. 

It follows from (|19|l and (|22|l that U in) is summable. Thus, boundedness of i?(n) 
follows from this and l|23|l (cf., e.g., 1211 Lemma 3.5]). This concludes the proof. □ 

The theorem above implies the stability of the number a and of the sets AiN), 
i3i(T), A under suitable power-decaying perturbations of the potential. On the 
other hand, the measure of the sets /i(i?i(T)), piA) and the Borel transform Fiz) 
may change after such a perturbation. In particular, it is possible that piA) ~ 
for the perturbed operator in part (b) of Theorem ^ Thus, the bounds (|10|l and 
H12|l are in general not stable. Of course, we still get a dynamical bound for the 
perturbed model. For example, we have the following consequence of Theorem^ 
and Theorem 12 

Corollary 3.1. Assume that for some energy Eq and some constant Ci, the trans- 
fer matrices T associated with A + y satisfy ||T(7i, tti; i?o)|| < Ci|7i — 77i|" for every 
n,m € Tj with nm > 0. Assume further that, for some e > 0, the perturbation W 
satisfies \Win)\ < C2\n\~^~'^°'~'^ for every n G Z. Then we have for the operator 
A + V + W, 

I3I ip) > 
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for every p > 0. 

Proof. By Theorem |5J we have that the transfer matrices T' associated with A + 
V + W satisfy \\T'{n, m; £^o)|| < C\n — rn\" for every n,m ^ Z with nm > 0. Then, 
an inspection of the proof of Thcorcmnshows that this suffices to prove the bound 
Hll|l which yields 

< \X\P^ > (T) > CT'Slr-i 

and the assertion of the corollary follows. More precisely, one can work indepen- 
dently on the two half-lines and hence needs bounds on \\T'{n, m; £')|| only for the 
case where n, m have the same sign. □ 

4. A Class of Pseudo-Sparse Potentials 

In this section, we study a class of "sparse" potentials which includes various 
substitution models and the prime model. These potentials are not all sparse in the 
standard sense, but the point is that the class we discuss contains sparse potentials, 
and also a number of other potentials that have been considered before and which 
can be studied within the same framework. 

Let us consider the case where the potential V is defined on the half-line N and 
takes on two values a, 6 e M. We assume the following for n large enough, that is, 
for n > N: 

(51) Occurrences of b are always isolated, that is, if V{n) = b for some n, then 
V{n- 1) = V{n + 1) = a. 

(52) The value a always occurs with odd multiplicity, that is, if V{n) = V{n + 
k + 1) — b and V{n + j) — a, 1 < j < k, then k is odd. 

Sparseness in this context refers to the 6's being isolated and the results below 
holding for arbitrarily long gaps between consecutive 6's. However, some of the 
concrete applications — for example the applications to substitution models — will 
not be sparse in a traditional sense. 

We can prove the following: 

Theorem 3. Suppose V : N {a,b} C M is a potential satisfying (SI) and (S2) 
above. We have for every p > 0, 

Proof. Up to an initial piece, the transfer matrices are given by products of matrices 
of the following form: 

T{a,Ef^+^ and T{b,E). 

Let Eo ~ a. Then 

T{a,Eo)^'+^ = {T{a,E^f)'T{a,Eo) = {-if T{a, Eo) ^ ±T{a, Eo). 
Up to sign, this gives rise to powers of 

Clearly, such powers satisfy a bound which is linear in the number of factors. Thus, 
the claim follows from Ullf) . □ 
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Remark. We can apply Corollary 18.11 and obtain that the dynamical bound in 
TheoremOlis stable with respect to perturbations W obeying |Vl^(n)| < C2n^'^~'^ 
for some fixed £ > and every n € N. Similarly, we have stability with respect to 
power-decaying perturbations for all the dynamical bounds that will be shown in 
this section and we will not make this explicit for each one of them. 

Let us now discuss the case where the a's occur with even multiplicities. That 
is, we assume for n large enough, 

(S3) The value a always occurs with even multiplicity, that is, if V{n) = V{n + 
k + I) ^ b and V{n + j) ~ a, I < j < k, then k is even. 

In this case we can prove a dynamical bound even without assuming the sparse- 
ness condition (SI). However, we need that \a — b\ is not too large. Namely, we 
have the following result: 

Theorem 4. Suppose V : N —i- {a, 6} C K zs a potential satisfying (S3) above. 

(a) If \a — b\ < 2, then for every p > 0, 

(b) // |a — 6| = 2, then for every p > 0, 

Proof. The argument proceeds in a way similar to the proof above. Again, up to an 
initial piece, the transfer matrices are given by products of matrices of the following 
form: 

T{a,Ef^ and T{b,E). 

Again, let £'o = a. Then 
and hence 

T{a,E^f' = {T{a,Eo)'y = {-!)'= ±1. 
On the other hand, r(6,i?o) is elliptic when |a— 6| < 2 and parabolic when |a — 5| = 
2. Thus, in the former case, products of matrices of the form T(a, i?)^' or T{b, E) 
remain bounded, while in the latter case such products satisfy a bound which is 
linear in the number of factors. The claim thus follows from □ 

Let us note that a result like part (a) of Theorem 0| is implicitly contained in 
|17| . where mainly random polymer models are studied. 

It is clear that whole-line analogs of the above theorems hold. In this case, we 
need (SI) and (S2) or (S3) to hold for \n\ large enough. 

More importantly, these results cover a variety of seemingly very different cases: 
First consider the period doubling Hamiltonian, which was already discussed in 
On the alphabet A = {a, 6} C R, consider the period doubling substitution 
S{a) = ab, S{b) = aa. Iterating on a, we obtain a one-sided sequence 

u — abaaabababaaabaaab . . . 

which is invariant under the substitution process. Define the associated subshift ilpd 
to be the set of all sequences over A which have all their finite subwords occurring 
in u. Here, we can consider either one-sided or two-sided sequences. This does not 
matter for the results in this paper, but we remark that for substitution models, 
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one generally considers the two-sided case. For lo G ripd, we define the potential 
by Vujin) — ij-'ri- It is easy to check that each satisfies (SI) and (S2) (even 
for every n € Z) and hence an application of Theorem |31 allows us to recover jSj 
Theorem 3]. However, we can prove a more general result. Consider, for example, 
substitutions of the form 

(24) S(a) = c?^-^h, S{b) = k,l>l. 

The case fc = 1, ^ = 1 corresponds to the period doubling case. The potentials 
generated by a substitution of the form (Pl|l (by generating a one-sided fixed point 
and passing to the associated subshift, as in the period doubling case above) are 
easily seen to obey (SI) and (S2). On the other hand, substitutions of the form 

(25) S{a) = a'^'^b, S{b) = a^', k,l>l 

give rise to potentials satisfying (S3) and hence Theorem 0] applies in these cases. 
Thus we may state the following: 

Corollary 4.1. (a) Let S be a substitution of the form (|24|l . f2 the associated 
subshift, and for lu Q, let V^(rt) — ujn, n e Z. Then, for every lo ^ U,, the 
potential gives rise to an operator satisfying 

{p) > ■1'°'^ ^^^'^y p> 0- 

(b) Let S be a substitution of the form (|25|1 . f2 the associated subshift, and for 
G ri, let Vcj{n) — LUn, 71 G Z. Then, for every lu € fl, the potential gives rise 
to an operator satisfying 

Psi (-P) ^ -P ^ 1 /""^ every p > if \a — b\ < 2 

and 

Ps, {p) > for every p > if \a - b\ ^ 2. 

Consider the following class of substitutions: 

(26) Sia) = a"'b'', S{b) = a. 

The case m — n — I gives rise to the Fibonacci substitution. Hence, the substi- 
tutions in (|2(j|l are usually called generalized Fibonacci substitutions. If n = 1, 
the resulting potentials are Sturmian and will be discussed in this more general 
context in a later section. Here, we restrict our attention to the case n > 2. These 
substitutions and the associated Schrodinger operators were studied, for example, 

in [22132 EH. 

If n is even, it is easily seen that each Vu, satisfies (S3) with the roles of a and b 
interchanged, that is, &'s always occur with even multiplicity. Thus, we can derive 
a dynamical bound for the associated operators by applying Theorem 0] 

If n is odd, the model satisfies neither (S2) nor (S3) but we can nevertheless 
employ a similar argument. As a warmup, let us consider the case n = 3 (the 
special case m = 1, n = 3 is usually called the nickel mean substitution). Then the 
transfer matrices are given by products of matrices of the following form: 

T{a,E) and T{b,Ef. 



QUANTUM DYNAMICS IN ONE DIMENSION II. 



11 



Let Eo^b+1. Then 
and hence 

T{b,Eof = -/. 

This would allow us to prove bounds on (p) in the same way as in the proof of 
Theorem 0] 

Let us now turn to the case of a general odd n > 3. Here we can extend the 
above idea and prove a result which applies to the substitutions in (|26|l with n odd 
but which is much more general. Denote 

(S4) There is some odd k > 3 such that the value b always occurs with a mul- 
tiplicity which is a multiple of k, that is, if V{n) = V{n + I + 1) = a and 
+ j) = I < j < I, then / = mk for some m e N, 

Then we can prove the following: 

Theorem 5. Suppose V : N ~f {a,b} <Z R is a potential satisfying (S4). Then 
there is a set £ dM. of cardinality k — 1 such that for every E Cz £, we have 

(a) // |a — ii'l < 2, then for every p > 0, 

(b) // \a — E\ ~ 2, then for every p > 0, 

Proof. In this case, the transfer matrices are given by products of matrices of the 
following form: 

T{a,E) and T{b,Ef. 
It suffices to exhibit fc — 1 energies £^o with 
(27) T{b,EQf^±I. 

This can be seen as follows: The matrix T(5, E^ is the monodromy matrix of the 
constant potential V{n) = b, regarded as a fc-periodic potential. This gives rise 
to an operator with fc — 1 gaps. However, since the operator with this potential 
has spectrum [6 — 2,6 + 2], all these gaps are degenerate. Every degenerate gap 
corresponds to an energy where the monodromy matrix is equal to ±/, hence there 
are exactly fc — 1 energies Eq for which we have H27|l . □ 

Putting everything together, we obtain the following result for the models gen- 
erated by substitutions from (|26|l : 

Corollary 4.2. Let S be a substitution of the form (|26|) . fl and the V^j 's as above. 

(a) Ifn>2 is even, then for every w G il, the potential V^j gives rise to an operator 
satisfying 

{p) > P — 1 for every p > i/ |a — 6| < 2 

and 

(P) ^ ^-5— ^""^^y P > if \a- b\ = 2. 
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(b) If n > 3 odd, then T(b,Eoy^ ~ ±/ has n — 1 solutions Eq Cz M. and for each 
such solution Eq, we have that for every cj G fi, the potential gives rise to an 
operator satisfying 

(p) > P - 1 for every p> if \a - Eq\ <2 

and 

Ps, (p) > f""^ ^""^^y P > «/ |a - E'ol = 2. 

The final substitution model we consider is the following: 



(28) S{a) = a'^b"', S{b) = 6"a'". 

The case m = n = 1 gives rise to the Thue-Morse substitution. Hence, the substi- 
tutions in (|28|l are usually called generalized Thue-Morse substitutions. They were 
considered, for example, in [3^. If at least one of m,n is even, (S3) holds and we 
can apply Theorem^ In the remaining case, where both m and n are odd (and at 
least one is > 3), (S4) holds and we can apply Theorem^ Thus, for models gener- 
ated by generalized Thue-Morse substitutions, we obtain the following dynamical 
bounds: 

Corollary 4.3. Let S be a substitution of the form (|28|) . Q and the 's as above. 

(a) // at least one of m, n is even, then for every w G the potential VL, gives rise 
to an operator satisfying 

Pg^^ ip) > P — 1 for every p > if \a — b\ < 2 
and ^ 

Psi (P) ^ for every p > if \a ~ b\ ^ 2. 

(b) // we have m > 3 odd, then Tib, i?o)™ — has to — 1 solutions i?o € K and 
for each such solution Eq, we have that for every uj € Q, the potential gives rise 
to an operator satisfying 

Psi (p) > P - 1 for every p > if \b - Eo\ < 2 
and ^ 

Ps, (p) > for every p > if \b - Eo\ = 2. 

An analogous result holds if we have n > 3 odd. 

(c) If m = n = \, then 

/^a^ (p) > P ^ 1 for every p > 0. 

Part (c) was proved in [S] and is stated for completeness. One might expect 
the bound (p) > p — 1 to hold always. In fact, the paper claims, for every 
choice of m, n, a, b, the existence of an energy, where the transfer matrices remain 
bounded. However, the argument given in that paper is incomplete and it would 
be interesting to prove or disprove this claim. 

Next, we consider the prime Schrodinger operator -ffprimo on ^^(N) whose poten- 
tial is given by 

a if n is not prime. 



Vprimo(n) - I ^ ^ p^j^g_ 

This operator was studied, for example, in |1(JI I3U| . Based on numerics and heuris- 
tics contained in these two papers, one may expect the following: On the one hand, 
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for almost every energy there is an solution to i/primo'/' = Eip, that is, when 
one varies the boundary condition at the origin, one gets pure point spectrum for al- 
most every boundary condition. On the other hand, the model displays non-trivial 
transport for every boundary condition. We will confirm the latter below (the proof 
discusses only the case of a Dirichlet boundary condition, but it readily extends to 
every other boundary condition). Let us briefly discuss the flrst point. It is natural 
to view Vpiimc as a sparse potential. In fact, this point of view was proposed in 
|1U| . However, the current methods in the spectral analysis of models with sparse 
potentials (see, in particular, |2l)[ \2<i\ ) are clearly insufficient to conclude anything 
for the prime model. We regard this as an interesting problem and refer the reader 
also to 29 for further motivation to consider models of moderate sparseness. 

Let us now turn to a dynamical result for the prime model. Clearly, (SI) and 
(S2) are satisfied for n large enough. Hence we get: 

Corollary 4.4. For every a, 6 G M, the operator i/piimc satisfies 

(p) ^ ^-y- /o'' P > 0. 

Finally, we discuss a model which is sparse in the standard sense. Namely, pick 
some integer 7 > 2 and define Uk = 7*^ for fc S N. Let K!parse(") = b if n = Uk for 
some k and Kiparso('T') = a otherwise. Schrodinger operators with potentials of this 
kind were studied in |35| . Clearly, when 7 is even, all n^'s are even, and when 7 is 
odd, all Tifc's are odd, so we have (SI) and (S2). Thus, Theorem O applies and we 
get 

Corollary 4.5. For every a, 6 G R and 7 G N \ {1}, the potential T4paiso gives rise 
to an operator satisfying 

{p) > ■1'°'^ ^'"^'^y p > 0. 

This can be improved if 7 ^ |a — 6|: 
Proposition 4.6. Let 

^ 2 log ^2+ (a -6)^ 
log 7 

Then the potential V^parsc gives rise to an operator satisfying 

_ p — 1 — Aiy 
Ps, (p) > — — for every p>0. 

J- 4-7/ 



Proof Write C(a, b) = ^2 + (a -5)2. Then 

\\T{a,E = af^+^T{b,E = a)\ 

For dnjrn — < k < n : V{k) = 6}, we have dn^m < log \n — m|/ log7 and hence 

\\T{n,m;E = a)\\ < C{a,bf"-"^ < C(a, 6)'°^ l«-™|/iog7 ^ |n - m|'°g^(''''')/'°s''. 



-1 
a ^ b 



< C{a,b). 



This yields the assertion. 



□ 
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5. A Hierarchical Model 
The hierarchical model is defined through the potential 

(29) V{n) ^ Xfiordn), 

where / is some real function and ordn is the number of factors 2 in the prime 
decomposition of n. The sequence ()29f) has some nice symmetries. Because 
ord (— n) = ordn for all n and ord {I ■ 2™ + k) = oidk for m > 1, alH and |A;| < 2™, 
analogous identities hold for V. In particular, 

(30) V{1 • 2™ + fc) = V{k) ^ V{-k) = V{1' • 2™ - fc) 

for any I and I', m > 1 and |fc| < 2™. The Schrodinger operator with such a 
potential appeared first in the works 26 and with the special choice 

m— 1 

f{m) = 

where i? is a positive constant. The advantage of this choice is that in this case, 

x,n = tiM„^{E) = ti-T{2"\0;E) 
satisfies an autonomous difference equation |31j . 

(31) Xm+l = - '2 + RXmiXm - + 2) , m>l. 

The above recurrence and the symmetries (|3U|I made it possible to obtain many 
rigorous results about the spectrum of the corresponding Schrodinger operator. A 
detailed mathematical study of this model was carried out by Kunz et al. |23j . 
Among other things, it was shown that for every i? > 0, the spectrum is a Cantor 
set, and for i? > 1, it is purely singular continuous. From the point of view of the 
present article, it is interesting that a countable infinite set of exceptional energies 
in the spectrum could be identified explicitly. The 2™ zeros Emk, 1 < fc < 2™, 
of Xm{E) are simple and Xm = implies Xm+i = — 2 and x^+i — 2 for / > 1; 
compare H31(l . From this it was possible to show that E^k, for to > and 1 < k < 
2™, are lower (resp., upper) gap-edges in the spectrum of if A > (resp., A < 0) 
and they are dense in the spectrum. For the corresponding gap-edge states, the 
following result was obtained (Proposition 15 in |23p. 

Proposition 5.1. Let Xm{E) = and let ip be a solution of Hip = Eip. 

(i) If ifjifi) = 0, then ip(k + 2™+'^) = —ipik) for every integer k. 

(ii) If ip{Q) ^ 0, then i(j{2l ■ 2™) = (— 1)'?/>(0) and asymptotically, as I ~* oo, 

(32) V((2/ + 1)2") - ^(2") X {-iy+'X.nmfRil) 
where 

{2^' , R<2 
Mog2^ R = 2 

( 2_Yi ;logo -R. p 9 

\r) 2(fl-l)(fl-2)' ' ' -n- ^ ^• 

Here Xm — XR"^ Xm-i{E) ■ ■ ■ xo{E) , ei G [0,1) is the fractional part of log2 I 
and X means equality in the leading order of I . 
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We use this proposition to prove the following theorem. 
Theorem 6. For every A ^ and R > 0, 

where 

a = a{R) — max{l, log2 R}. 

Proof. We apply Proposition 15.11 with m = for which it provides the precise 
asymptotic form of the solutions. Because xo{E) = E, these belong to = 0. Let 
ipD and ip-^ be the two solutions defined by the initial values 

(34) Vd(0) = Vn(1) = 0, V^d(1) = ^n(0) = 1. 

According to part (i) of Proposition 15.11 V'd is a periodic solution with period 4, 
namely 

(35) Vd(20 = 0, i/;d(2/+1)-(-1)'. 
On the other hand, 

(36) V^n(20 - (-1)', ^n{2l + l):<{-iy+^XfR{l). 

Equations l|35|) and (|36|1 permit us to compute the asymptotic form of T{n,m;0). 
Because of V{—n) = V{n), it sufhces to consider n > m > 0. In what follows, we 
use the simplified notation T{n,m). Let ^'(ti) — {'ip^{n + 1) 'ip^{n))'^ for « = 0, 1. 
Then r(n, 0) = (^'d('^) ^n(?^))- The determinant of any transfer matrix being 
unity, the inverse is easy to compute. We find 

(37) T(n,m) = r(n,0)T(m,0)-i 

ijTtin + l) Vn(?^+1)W i'^im) --0N(m + l) 
V'D(n-) V'n(") / V -V'D(n^) VD(?n + l) 

With the short-hand notation 

i^(0 = (-l)VN(2/ + l), 
equations H35|l . (|36|l . and H38|l then yield 

1 F{l)~F{k) 
1 

-1 

1 F{l)-F{k) 

F{l)-F{k) -1 
1 

-1 
F{l)-F{k) -1 

All these matrices have the same norm. Denoting the Hilbert-Schmidt norm by 
II • II2, for n = 2Z, 2Z + 1 and m = 2fc, 2/c + 1, we have 



(38) 





T(2;,2fc) 




-1)'=+' 




T(2« + 1,2/c) 


= (- 


-1)'=+' 




r(2/,2A:+ 1) 


= (- 


-|^-jfe+;+i 


(39) 


T(2; + l,2fc + l) 


= (- 


^■jfc+'+i 



||T(n,m)|| < ||T(n,m)||2 = ^2 + [F{1) - F{k)Y ^ ^2 + X^/rH) - fnik)]^ 
Therefore, 

||T(n,m;0)|| <2A/fl(n/2) 
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for any n large enough and m < n. If i? 7^ 2, the assertion of the theorem obviously 
follows from the definition H3^^|l of /r and Theorem ^ If i?, = 2, we note that for 
any e > 0, 

||r(n,m;0)|| < Ani+' 
if n is large enough. Therefore, by Theorem ^ 

for any e > and, thus, for e = as well. □ 

Remark. The proof shows that we can apply Corollary 13.11 and obtain that the 
dynamical bound in Theorem is stable with respect to perturbations W obeying 
< C2\n\^^^^°'^'^ for some fixed e > and every n e Z. 

We note that instead of to = 0, we could have used Proposition 15.11 with any 
TO > and any zero of Xm{E). This holds because of the following: 

Theorem 7. For any X ^ 0, R > 0, m > 0, and k E {1, 2, . . . , 2™}, there exists a 
positive number C\,ii{m, Emk) such that for any n> n' > 0, 

\\T{n,n';E^k)\\ < Cx,R{m, E^k)fR{2-"'-^n) . 

Proof. We fix TO > and a zero Emk of Xm. From equation 1)38(1 it is clear that 
we have to bound the two particular solutions ((34|l of Hi/; = Emkip. According to 
Proposition 15. II V'd is 2'"+^-antiperiodic and, thus, bounded. On the other hand, 

(40) ^n(2/-2™) = (-1)', V'N((2/ + l)2'")-VN(2™)x(-l)'+iA„Jfl,(0. 
Thus, the task is to bound ip't^{n) in the intervals 

(41) 2/ • 2™ < 7i< (2/ + 1)2™ and {21 + 1)2" < n < 2{l + 1)2" . 

To proceed with the proof, let us recall equation (3.29) of 23 , according to which 

■00(2") = Xm-l ■■■Xq 

for any energy. Thus, -00(2™) 7^ in the present case {E — Emk), for otherwise 
Xi = for some i < m would imply \xj\ =^ 2 for every j > i, contradicting Xm — 0. 
Then uo :— ?/'d/V'd(2'") is a solution of the Schrodinger equation satisfying the 
boundary conditions uo(0) = 0, ^0(2™) = 1 and, according to Proposition 15.11 
wo(fc+2™+^) = — Mo(fc) for any k. From the general theory of second-order difference 
(differential) equations, it follows that there exists a linearly independent solution 
Ml with boundary values ui(l) = 1, ui{2"^) = and that we can write ■i/'N for 
<n < 2'" in the form 

V'N(n) = i/;N(2")Mo(n) + ^N(0)Mi(n). 

Next, we observe that ui can be expressed in terms of uq. Indeed, from equation 
(I3U|I we can see that the sequence V{1), . . . , V{2"^ — 1) is a palindrome, 

l/(2™"i - fc) = T/(2™-i + fc) fc = 1, . . . , 2"-i - 1 

and, hence, 

ui(n) =Mo(2"-n), n = 1, . . . , 2" - 1. 
Furthermore, the translational symmetry of the potential, 

{v{i • 2™ + 1), . . . , viii + 1)2™ - 1)) = (y(i), . . . , y(2™ - 1)), 



QUANTUM DYNAMICS IN ONE DIMENSION II. 



17 



valid for any I, implies that the translates of uq and ui can be used to give i/'N in 
each of the intervals H41(l . Altogether we find 

ipNin) = ?An((2/ + l)2™)uo(n - 21 ■ 2™) + ?/'n(2/ • 2™)uo((2Z + 1)2™ - n) 

if 2/ • 2™ < n < {21 + 1)2" and 

V;n(7i) = VN(2a + l)2'")uo(« - (2; + 1)2"^) + ?/'n((2Z + 1)2™)mo(2(/ + 1)2"* - n) 

if (2Z + 1)2" <n< 2{l + 1)2". Together with in both intervals, 

\^n{n)\ < j^iM(|^^((2/ + i)2")| + 1). 
Since I < n/2™+^, we obtain that for n large enough 

I^nWI < |^iM(|A,„|/^X2""-^n) + |V'n(2™)| + 1). 
Due to H38|l . the assertion of the theorem follows from this bound. □ 

6. Sturmian Potentials 

In this section, we discuss dynamical bounds for the standard one-dimensional 
quasicrystal model which is given by a Schrodinger operator on the whole line whose 
potential is given by 

(42) y(n) = Xv^^g{n), where ve{n) = X[i-ujA){nuj + 9 mod 1), 

where A 7^ is the coupling constant, lu e (0, 1) irrational is the rotation number, 
and 9 G [0, 1) arbitrary is the phase. For more information on this family of 
operators, we refer the reader to the survey articles 4, 33 . 

It is well known, and easy to see, that the spectrum of the operator Hx ^^ g with 
potential V from H42I) is independent of 6', that is, for every A, lo, there is a set Y.\_^ 
with o-(i?A,a;,e) = '^\,u> for every 9. 

Consider the continued fraction expansion of w, 

1 



1 

ai H 



1 

02 H 

as H 

with uniquely determined an G N (cf. |18)). The associated rational approximants 
Pk/qk are defined by 

Po = 0, _pi = 1, Pk ^ a-kPk-i + Pk-2, 
qo = 1, gi = ai, qk = CLkqk-i + qk-2- 
The number tu is said to have bounded density if 

1 " 

(43) d{u!) = limsup — Ofe < 00. 

n — ►oo , 

k—1 

The set of bounded density numbers is uncountable but has Lebesgue measure zero. 
The following was shown in (see also ^1] for the case of zero phase): 



18 



D. DAMANIK, A. SUTO, S. TCHEREMCHANTSEV 



Theorem 8. Suppose uj is a bounded density number. For every \, there is a 
constant C such that for every 9, every E G Sa,^, and every n,m G Z, we have 

(44) \\Tx.,u.An, m- E)\\ < C\n - m\"^^^'^\ 
with 

(45) a{X,uj) = D ■ d{uj) -logCx, 
where D is some universal constant, C\ is given by 



(46) Cx = 2+ v/8 + A2, 

and d{uj) is as in (|43() . 

This yields the foUowing: 

Corollary 6.1. Let uj be a bounded density number. Then, for every X,9, the 
operator H\^^^g satisfies 

n- r \ P " 3a(A, Lij) . 

Ps. (P) > — — 7T — ^ for every p>0, 

with a{X,Lu) given by (|45|l . 

Since /x(S]a,w) = 1, this is an immediate consequence of (I10() . This bound is 
better than the corresponding result in [S] (which follows from bounding from 
below \Bi{T)\). One should stress that as opposed to all the other examples dis- 
cussed earlier, the dynamical bound in Corollary 16. II is not stable with respect to 
perturbations of the potential. This is due to the fact that fi{Yi\^^) may vanish for 
the perturbed measure. However, by Corollarv l3.1l we have the following result: 

Corollary 6.2. Let uj be a bounded density number and let A be arbitrary. If 
a(A, uj) is given by (|45(l and W satisfies 

\W{n)\ < C2(l + |n|)-i-2a(A.c^)-e j-^^ g^g^^ ^ g ^ 
for some £ > 0, then, for every 6, the operator H\ i^^g + W satisfies 

. , p — 1 — 4q!(A, w) 
Pa, \P) > — 7-, — 7T — ^ — for every p> 0. 
^ 1 + a(X, ij) 

As in the case uj = (\/5 — l)/2 and 9 — 0, studied in [H|, it is possible to 
improve this lower bound somewhat by exhibiting a suitable set A{N) (stable under 
perturbation), studying its Lebesgue measure, and applying (j^. The set A{N) will 
again be given by the spectra of suitable periodic approximants, and the Lebesgue 
measure can again be bounded through a fine analysis of the trace map, akin to 
what is done in |51 ll9[l77| : compare also |^. We leave the details to the interested 
reader. 
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